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Abstract. Production and decay angular distributions were extracted from measurements of exclusive 
electroproduction of the ^"(770) meson over a range in the virtual photon negative four-momentum squared 

0.5 < Q'^ < 4 GcV^ and the photon-nuclcon invariant mass range 3.8 <W< 6.5 GcV. The experiment was 
performed with the Hermes spectrometer, using a longitudinally polarized positron beam and a "^He gas 
target internal to the HERA e* storage ring. The event sample combines mesons produced incoherently 
off individual nucleons and coherently off the nucleus as a whole. The distributions in one production angle 
and two angles describing the 7r"'"7r~ decay yielded measurements of eight elements of the spin-density 
matrix, including one that had not been measured before. The results arc consistent with the dominance of 
helicity conserving amplitudes and natural parity exchange. The improved precision achieved ai A <W <7 
GcV, in combination with other data at > 7 GcV, reveals evidence for an energy dependence in the 
ratio R of the longitudinal to transverse cross sections at constant . 
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1 Introduction 



The study of exclusive photoproduction and leptoproduc- 
tion of vector mesons (leaving the target nucleon(s) in- 
tact) has provided information both on the hadronic com- 
ponents of the photon, and on the nature of diffractive 
scattering. Theoretical models typically describe the ex- 
clusive production of light vector mesons as occuring via 
the fluctuation of the real or virtual photon into a quark- 
antiquark pair (or off-shell vector meson), which is scat- 
tered onto the mass shell by a diffractive interaction with 
the target. The corresponding tree-level diagram is shown 
in Fig. |], and a list of kinematic variables describing the 
reaction is given in Table |l]. The positron emits the virtual 
photon with energy v and four- momentum squared — Q^, 



which then scatters from the target by the exchange of 
particles in the ^-channel. 

The decay-angle distribution is influenced by the trans- 
fer of the virtual photon's spin to that of the vector me- 
son and can therefore provide a stringent test of the re- 
action mechanism. Electroproduction data are now avail- 
able with modern precision in the range of the photon- 
nucleon ccntrc-of-mass energy W reaching from low {W < 
4 GeV) |,|,|,|, through moderate {6 < W < 20 GeV) i, 
|fy|l , and recently to very high (50 < < 140 GeV) |, 
lOrri) values, and for values of up to 30GeV^. A vari- 
ety of competing models mentioned below have achieved 
some success in reproducing the observed dependences of 
the virtual photoproduction cross section on Q^, W, and 
t. The polarization data presented here cover a range in 
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Table 1. Kinematic variables used in the definition of exclusive diffractive vector meson events. (Definitions involving energies 
or three- vectors are given for the laboratory frame.) 



Variable Description 



k = 


{E,-k) 
{E',k') 




4-momentum of incident positron 


k' = 




4-momentum of scattered positron 


P = 


(Mn,0) 




4-momentum of target nucleon N (ignoring Fermi motion and binding) 


q = 


{v, q) = k — k' 




4-momentum of virtual photon 








Photon virtuality 




e 




Positron scattering and azimuthal angles in the lab 






vlEf\ 


((^ic > defines upper half of spectrometer) 


£ = 


{X-vlE)l\\-vlE^\{ 


Virtual photon polarization parameter (valid for <C Q'^ -C E^) 


w = 


--\'y1 + P)\ = (A^N + 2Mni^ - 


Total photon-nucleon centre-of-mass energy 


X = 






Bjorken scaling variable 


y = 


vjE 




Fractional leptonic energy transfer 


Ph± 






3-momentum of positive or negative hadron 




= y/m| + p 2^ 




Hadron energy assuming it is a pion 


pv = 


= (£v,pv) = 


-,Ph+ +Ph-) 


Candidate fp 4-momentum 








Candidate fp mass 


My 


= \/(^' + 9-Pv)' 




Invariant mass of the recoiling baryonic system 


t = 


(g-pv)'' <0 




Squared 4-momentum transfer to target 


to 






Minimum of \t\ for fixed Q^, M^^, and My 


t' = 


t - io < 




4-momentum transfer squared relative to the case 
where the meson and photon are coUinear 


AE 


= (M^ -M^)/2Mn = 


- £v + t/2MN 


Measure of energy absorbed by recoil system 



W that fills the gap between 4 and 7 GeV containing little 
previous data [|2|, and should therefore serve as a useful 
additional constraint on some of these models. 

The Vector Meson Dominance (VMD) model Q has 
traditionally been used to describe the reaction. It treats 
a photon with of order a few GeV^ or less as a superpo- 
sition of the lightest vector mesons that diffractively scat- 
ter from the target. As increases, the hadronic content 
of the photon as a quark-antiquark pair is resolved. More 
fundamental calculations including an explicit QCD-based 
treatment of the dynamics of the qq pair provide a unified 
treatment of the entire range for either the total 

or the longitudinal [jis) cross section. 



/ (^') 



y (^) 



(0 



.+ 



N(P) 



Y(P') 



Fig. 1. Tree level diagram of diffractive vector meson electro- 
production viewed as two-body virtual photoproduction. Four- 
momentum labels are indicated in parenthesis. 



In both VMD g and QCD-based models |T|,|T|,0: 
the phenomenological Regge theory is used to describe the 
diffractive reaction. At low energy {W < 4 GeV) the pro- 
cess is described in terms of the exchange of mesons (p, 
uj, /2. ap, Lu^, Pi ■,■■■), all lying on the same Regge trajec- 
tory h§. This so-called Reg geon exchange describes the 
decrease with W of the photoproduction cross section for 
real and quasi-real photons. Above W ~ 10 GeV, how- 
ever, the cross section increases slowly, and Regge models 
of the reaction involve the i-channel exchange of a color- 
less spinless object [^,16,1^, which has been identified 
with the soft Pomeron introduced to explain elastic pp 
scattering Q. 

The diffractive interaction has also been treated in 
perturbative QCD models of longitudinal production, 
which illuminate how the Reggeons arise at a fundamen- 
tal level. At the high x values corresponding to low W 
for moderate Q^, the diffractive interactions arise as the 
result of the exchange of two quarks Jl^ , corresponding 
to meson exchange in Regge theory. At low x (or high 
W), the interaction is assumed to involve the exchange of 
two gluons which has led to attempts to understand 
the phenomenological Pomeron in terms of nonperturba- 
tive exchange of gluons ||2^. At the highest Q^, the di- 
vergence of the nucleon's gluon content at low x leads to 
a strong W dependence ||l^,|l^,^,|2^. At the intermedi- 
ate energies (3.8 <W < 6.5 GeV) and photon virtualities 
(0.5 < < 4 GeV^) available to Hermes, both the quark 
exchange and the gluon exchange (or, in the Regge lan- 
guage, both Reggeon and Pomeron exchange) processes 
are expected to contribute. 
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For unpolarized targets and polarized beams, the dis- 
tribution in the production and decay angles provides in- 
formation on up to 26 independent elements of the spin 
density matrix ||2^, providing much tighter constraints 
on theoretical models than the total cross section alone. 
Furthermore, the angular distribution yields important 
model-independent insight into the reaction. For example, 
experiments have indicated that the helicity of the pho- 
ton in the 7*iV centre-of-mass system is approximately 
retained by the vector meson, a phenomenon known as 
s-channel helicity conservation (SCHC) Given the ap- 
proximate validity of SCHC, the angular distribution can 
be used to calculate the ratio i? = ctl/ot of longitudinal 
to transverse virtual-photon cross sections for exclusive 

meson production, without performing an experimen- 
tally difficult longitudinal-transverse (Rosenbluth) sepa- 
ration using two different beam energies. This ratio is ex- 
pected to be especially sensitive to details of the theoret- 
ical models. For example, in ref. it is asserted that 
R{Q^) is very sensitive to nonperturbative effects in the 
longitudinal quark-photon coupling. Recent calculations 
based on off-forward parton distributions have been used 
to describe the longitudinal cross section for exclusive p*^- 
production [^|2^,^ . Some features of these models im- 
ply that in a complete QCD description higher-twist ef- 
fects have to be large, a fact which should be kept in mind 
when interpreting the present data, especially in view of 
the relatively low Q^-domain of Hermes. 

Phenomenological models are often simplified by the 
assumption of SCHC. This assumption can be examined in 
calculations of spin observables including a detailed treat- 
ment of the underlying quark-gluon degrees of freedom. 
Violations of SCHC could arise from gluon ladders and 
quark loops ||2^; moreover, recent models of the quark-Po- 
meron coupling that address issues associated wit h g auge- 
invariance lead to non-zero spin- flip amplitudes |p8| . 

This paper reports the results from the f 995 Hermes 
data for the angular distribution of p°{770) production 
and decay, from which eight matrix elements related to 
the (0° spin density matrix were extracted. Some small 
matrix elements that are sensitive to violations of SCHC 
are thereby better constrained in this energy region, in- 
cluding some that are accessible only using the longitu- 
dinal polarisation of the lepton beam. These data also 
provide the most precise determination of R in this in- 
termediate kinematic region containing only one previous 
data set . In combination with existing data at higher 
energy, these new data allow the VF-dependence of R to 
be investigated in the region above W — 4 GeV, excluding 
the lower region where there are indications that the char- 
acter of the reaction changes quickly. Below W = A GeV, 
the cross section increases rapidly with decreasing W, and 
there appears evidence of significant resonance contribu- 
tions as well as a large phase difference between longitudi- 
nal and transverse production that suggests the presence 
of non-diffractive processes jl) . 

In Section |[ the experimental conditions, event se- 
lection and treatment of backgrounds are described. The 
representation of the general decay-angular distribution in 



terms of a set of matrix elements is presented in Section 
^. The extraction of the matrix elements from the data 
is described in Section ^. The results are presented and 
discussed in Section ^, and the paper is summarized in 
Section ^. 



2 Experiment 

The experiment was performed using a gaseous ^He tar- 
get contained in an open-ended cell coaxial with the po- 
sitron beam at the HERA ep collider. During 1995 HERA 
produced a 27.5 GeV positron beam with an average lon- 
gitudinal polarization of 0.48 ± 0.03. Luminosities up to 
2x10'^^ nucleons-cm~^s~^ were achieved. Scattered posit- 
rons and coincident hadrons were detected in the Hermes 
spectrometer, of which a detailed description can be found 
elsewhere p9| . The spectrometer is equipped with front 
and back drift chambers before and after a 1.3 T-m dipole 
magnet. This tracking system provided a momentum res- 
olution Ap/p « 1-2% in the 1995 running period. An iron 
plate that shields the circulating beams from the bending- 
magnet field separates the spectrometer into identically 
instrumented upper and lower halves. The positron trig- 
ger is produced by a timing hodoscope in coincidence with 
an electromagnetic calorimeter system; this system com- 
prises a totally absorbing lead-glass array and a preshower 
detector consisting of a passive lead radiator and a scintil- 
lator hodoscope. This trigger required a deposit of at least 
4 GeV in the lead-glass array. A gas threshold Cerenkov 
detector and a transition radiation detector provided ad- 
ditional positron/hadron separation. 



2.1 Selection of Events 

On-line hardware monitoring systems ensured that the 
only events considered were those recorded when the appa- 
ratus was fully operational. Events with exactly one iden- 
tified positron and two oppositely charged hadrons were 
selected. Fiducial cuts were placed on the tracks to ensure 
that they all originated from the target and fell within the 
spectrometer acceptance. Positrons and hadrons were dis- 
tinguished with a likelihood analysis that combined the re- 
sponse of the four particle-identification detectors (see [ p9| 
and references therein). 

The kinematics of the detected positron define the 
four-momentum q and polarization parameter e of the 
virtual photon. The kinematic constraints applied during 
analysis are 

W > 3.8 GeV, 
Q2 > 0.5GeV^ 
y < 0.85, 
I sin 6*0 sin (/)o I > 0.04. 

The invariant mass Af^^ of the hadron pair was recon- 
structed under the assumption that both hadrons were 
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Fig. 2. Mtttt distribution for exclusive diffractive events (filled 
circles) . The solid line is th e vi eld from the Monte Carlo simu- 
lation, discussed in Section |4.l| . The peak at low Mtttt is due to 
the K^K~ decay mode of exclusively produced <^ mesons (see 



section 2.2) 



pions (see Fig. || 
the requirement 



The ^"^(770) mesons were identified by 



0.63 < M^^ < 0.91 GeV. 

The diffractive cross section is forward peaked for small 
t' ~ -p±{p°) (see Table 0): 

Here the slope parameter b is related to the RMS radii rp 



and Tn of the meson and nucleon or nucleus |30 



(2) 



Hence diffractive events are generally characterized by 
small values of t' leaving the gently recoiling target rem- 
nant well separated from the forward-going vector meson. 
The acceptance of the Hermes spectrometer precludes de- 
tection of the recoiling target or target remnant; however, 
because of the small energy width of the HERA positron 
beam (aE ~ 25MeV) and the good momentum and an- 
gular resolutions of the Hermes spectrometer, the energy 
AE absorbed by the undetected target remnant can be 
reconstructed (see Table |l| and Fig. ||) and used to reject 
events in which the target nucleon(s) do not remain intact. 
Exclusive diffractive events were consequently selected by 
adding the kinematic requirements 



AE < 0.7 GeV, 
-t' < 0.5GeV2. 



(3) 



Both coherent scattering from ■^He and incoherent scatter- 
ing from nucleons are contained in the peak at AE « 0, 
though the peaking in t' is stronger for coherent scatter- 
ing because of the larger size of the '^He nucleus com- 
pared with the nucleon (Fig. ^). Both contributions were 
included in the present analysis. The above experimental 
requirements were applied, to the extent applicable, in all 
of the results reported here, resulting in 4800 accepted 
events. 
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Fig. 3. AE distribution for candidate diffractive p" production 
events (circles), showing the exclusive peak near zero. The his- 
togram shows the estimated contribution from non-diffractive 
processes (see section 2.2.1). 
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Fig. 4. f'-distribution for exclusive production events, with 
a fit (solid line) of ( p^ to the data. The contributions from 
coherent and incoherent scattering are shown with dashed and 
dotted lines, respectively. 



2.2 Corrections for Background 

Background arises from both diffractive and non-diffrac- 
tive processes. Inelastic fragmentation, i.e. hard scatter- 
ing on partons followed by hadronization, is the primary 
source of non-diffractive background. The diffractive back- 
grounds include exclusive non-resonant tt+tt^ production 
in which the pions are produced without forming a res- 
onance, double-diffractive dissociation in which the tar- 
get remnant breaks up during diffractive vector-meson 
production, and exclusive production of to mesons. The 
K+K~ decay mode of the (j) meson is reconstructed with 
M^TT in the region below 0.63 GeV when pion masses are 
assumed for the kaons, and is thus excluded from the se- 
lected data sample. The pair production (7 — s- e+e~) con- 
tamination of the di-hadron sample resulting from mis- 
identification of the leptons is negligible. 
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2.2.1 Non-diffractive Background 

Non-diffractive processes constitute the dominant source 
of background, for which a correction was applied using a 
Monte Carlo simulation. Background events were gener- 
ated according to the non-diffractive deep-inelastic cross 
section using LEPTO, and the fragmentation was treated 
using the LUND model, both described in ref. The 
detector response was simulated with a GEANT-based 
Monte Carlo code ||2^. Three-track eN eh+h'X events 
were selected from these Monte Carlo events and subjected 
to the same analysis criteria as the data, thereby produc- 
ing a sample of non-diffractive background events. The 
Monte Carlo estimation was normalized to the data in 
the inelastic region defined by AE > 4GeV (see Fig. 
The experimental yields for AE < 0.7 GeV were corrected 
by 6% at small to 15% at large Q^, using this sample 
of background events. 



2.2.2 Diffractive Background 
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Fig. 5. A parameterization of the target excitation spectrum 
of double-diffractive dissociation (arbitrary units) . 



by the contribution to the region at larger AE by the dom- 
inant decay branch uj — > tt^tt^tt^. In principle, the small 
non-resonant contribution should be included in any the- 
oretical treatment that is compared to the data. 



Measurements of double-diffractive pp scattering have sug- 
gested that the mass spectrum of the baryonic system {Y 
resulting from dissociation of the target has the form 



Mi 



da 



1.8GcV^-(MN+m„)2 
1 



if Af2 < 1.8GeV2 
otherwise. 



This function is shown in Fig. ^ expressed in terms of 
AE = (M| - M^)/2Mn. Only 15% of diffractive events 
with target dissociation described by such a model would 
pass the AE cut used in this analysis, which is a bene- 
fit of the good energy resolution compared to that avail- 
able in previous fixed-target leptoproduction experiments 
cited herein. Recent measurements of diffractive /o° pro- 
duction with target dissociation at centre-of-mass ener- 
gies from 60 GeV to 180 GeV give aj-p^pOy /cr-yp^pOp = 
0.65±0.17, with a diffractive slope parameter bno — 2.1± 
O.TGeV"^ Q. Assuming that the Hermes acceptance 
is unaffected by target dissociation, this result implies 
that the contamination of the exclusive data by double- 
diffractive events is less than (6 ± 2)%. Moreover, no dif- 
ference in the measured longitudinal fraction of p° mesons 
was seen between the two processes This is not sur- 
prising as target dissociation is not expected to affect the 
helicity structure of the 7*p° vertex. Hence no correction 
or uncertainty was applied for this source of background 
since this contribution to the uncertainty is expected to 
be small compared to others. 

The measured distribution includes interfering contri- 
butions from both diffractive production of Ci;(783) — > 
TT~^Tr~ mesons (branching ratio 2.2%) and non-resonant 
exclusive production of tt+tt^ pairs, neither of which are 
distinguished here from diffractive p'^ production . How- 
ever, the Lu production cross section is down by an order of 
magnitude, and its contribution tends to be cancelled by 
the enhancement of the background subtraction induced 



2.3 Radiative Effects 

Calculated radiative corrections to the lowest order di- 
agram (Fig. |l]) have approximately a 20% effect on the 
measured exclusive vector meson production cross section. 
Radiative effects on polarization phenomena are typically 
much smaller. External radiative effects in the spectrome- 
ter are entirely accounted for in the GEANT-based Monte 
Carlo simulation. Because of the negligible target thick- 
ness and windowless cell, external bremsstrahlung before 
scattering is negligible. No corrections were applied for in- 
ternal radiative effects; they are mitigated by the cut on 
AE (^) , which ensures that no hard photon is radiated in 
the events used to reconstruct the production and decay 
angles. A systematic uncertainty was assigned for internal 
radiation, as discussed below. 



3 Angular Distribution of the Production and 
Decay 

3.1 The s-Channel Helicity Frame 

Analyses of previous data on exclusive vector-meson pro- 
duction have shown that only the s-channel helicity frame 
leads to matrix elements that are not strongly i-depend- 
ent [^|3^. The s-channel helicity frame is the centre-of- 
mass system of the j*N system, with the quantization 
axis of the p'~' opposite to the direction of the recoiling 
target. The reconstructed kinematics are essentially indis- 
tinguishable for coherent and incoherent scattering, cor- 
responding to the use of the nuclear or nucleon mass. The 
azimuthal production angle ^ (see Fig. ^) is defined in this 
frame as the angle between the lepton scattering plane and 
the p'^ production plane containing the momentum trans- 
fer to the target. The decay-angle distribution can then 
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Fig. 6. The s-channel helicity frame. 



be expressed in the p° rest frame in terms of the polar 
and azimuthal decay angles 6 and (j) relative to the same 
quantization axis, with the azimuth defined relative to the 
vector meson production plane (see Fig. The angular 
distribution W(cos 0, cA, is then completely character- 
ized by the three angles 6, (p, and <P. The spin state of 
the p*' is reflected in the orbital angular momentum of the 
pion pair, giving rise to the (cos 9, (p) dependence. The ^ 
dependence of W arises because the virtual photon polar- 
ization depends on 



3.2 Decomposition of the Angular Distribution 

A complete derivation of the decay-angle distribution can 
be found in ref. . Briefly, the decay distribution can be 
related to the polarization density matrix pw of the vec- 
tor meson expressed in the helicity basis, using the Wigncr 
£)-functions: 



W(cos 6*, </>,<!>) 



_3_ 



E 



,0,~<^)pAA'(<^)^i'.o( 



»• (4) 



AA' = -1.0,1 



The <P dependence appears only in pxy , which can be de- 
composed into 9 independent contributions characterized 
by the matrix elements p"y : a = 0, 1, 2, 3 represent trans- 
verse photons: unpolarized, the two directions of linear 
polarization, and circular polarization, respectively. Pure 
longitudinal photons correspond to a = 4. The remain- 
ing matrices a = 5. . . 8 are attributable to interference of 
longitudinal with transverse photons. Expressions can be 
found in Appendix A of ref. p3| giving the various p"y 



in terms of combinations of the Wick amplitudes for 
a photon of helicity 7 to produce a meson of helicity A, 
summed over nucleon spins. 

At fixed beam energy, it is not possible to explicitly 
separate the contributions from longitudinal and trans- 
verse photons, and the angular distributions are described 
in terms of the following linear combinations of the p"y : 



^\X' — 



1 + eR 



1 + eR' 



a = (1,2,3), 



r° , = VR , a =(5,6,7,8). 

Here R is the ratio of longitudinal to transverse cross sec- 
tions and e is the polarization parameter of the virtual 
photon. Equation (0) becomes p3 

- V2 3?(r'j'^) sin 29 cos (/) ~ r^t^ sin^ 6* cos 2<?!) 
-e cos 2$ (jl^ sin^ 9 + r^,^ cos^ 9 - V2 ^{rlg ) sin 29 cos 

—rl_i sin'^ 9 cos 2(f>] 



-e sin 2<P (^\/2 5(r^o) sin 29 sin (f> + 3(r^_;^ ) sin^ 9 sin 20 
-I- \/2e{l + e) cos <P (^r^^ sin^ 9 + r^^ cos^ 9 

-\/2 3?(r5^o) sin 29 cos - r^_i sin^ 9 cos 20 



+ y^2e{l + e) sin^ x 

'V2^{r%) sin 20 sin + 3(r?_i) sin^ 61 sin 20 



-Ph 



(V2^{rl„) sin 261 sin ^(rf^i) sin^ 6lsin20) 
-\/2e{l ~ e) cos^ x 

'\/2 3(rJo) sin 26* sin0 + 3(r[_i) sin^ 9 sin 20 



+ \/2e{l - e) sin^ x 

(rf^ sin^ 9 + rgp cos^ 9 - V2 !R.{r%) sin 29 cos ( 



-rl_i sin^ 9 cos 



20)]} 



(5) 



where Pf, is the longitudinal polarization of the lepton 
beam, and the approximation ^ has been made. 

The statistical precision of the present data set re- 
quired that the analysis be limited to only one-dimensional 
angular distributions. 



3.2.1 The cos9 Distribution 

Integrating (|^) over and averaging over ^ gives 

W(cos 9) = l[l~ rSo^ + (3rS4 _ l)cos20] . (6) 
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This distribution depends on only one matrix element Tqq , 
which can be identified as the longitudinal fraction of p'^ 
mesons. yV(cos0) is not constrained by SCHC. 



3.2.2 The ?/' Distribution 

The polarization angle -0 = — ^ is the angle between 
the decay plane of the and the plane of polarization of 
the photon (i.e. the lepton scattering plane). If SCHC is 
valid, the azimuthal dependence is a function of tj) only, 
and integration of (||) over cos 9 gives 



COS 



2^] 



(7) 



which can be interpreted as an interference pattern aris- 
ing from the amplitudes for production by the two linear 
polarization states of the photon. A positive value for r\_i 
indicates that the meson and photon spins are aligned such 
that decay pions are preferentially emitted in the lepton 
scattering plane. 



4 Extraction of Matrix Elements 

To account for radiative and instrumental effects on the 
measurements, a Monte Carlo event generator for exclu- 
sive diffractive production was developed and used in 
conjunction with the GEANT detector simulation. The 
matrix elements were extracted as those values that best 
fitted the Monte Carlo distributions to the data. 



4.1 Monte Carlo Event Generator 

The diffractive p° generator produces eN eNp^ events 
according to the leptoproduction double-differential cross 
section 101 



da 



1 



dQ'^dE' AttEE' 



-MQ^E'MQ'), 



(10) 



where the virtual photon fiux is defined according to the 
Hand convention 13811 



3.2.3 The d> Distribution 



EriQ^E')^ 



1 E' 2 



47r2 2Mn Q'^ E l-e' 



(11) 



Integrating (||) over cos and averaging over gives 
W(0) = 

^ [l-2rOiiCOs20 + Pbx/r^3(rf_i)sin20] . (8) 



The azimuthal decay angle (j) is sensitive to the linear po- 
larization of the p*^ and (^ produces only terms of the 
form e'^^-^')"^ with |A - A'| = 2. This distribution mea- 
sures any SCHC-violating couplings of the linear polar- 
ization (A = 1,A' = —1) of the vector meson to trans- 
verse unpolarized [a = 0) and circularly polarized (a = 3), 
and longitudinal {a = 4) photons. If SCHC holds, then 
>V(0) = l/27r. 



Using VMD, the virtual photon production cross section 
can be expressed_in terms of the real photoproduction 
cross section as 



a(Q') = 



1 



(1 + QVMIY 



a(Q2=o), (12) 



a(Q2 = 0) = A 



2Afn 



(13) 



with = 29.4 pb • GeV and = 9.5 pb from a fit to 
data from ref . Q . 

The p*^ mass selection was weighted according to a 
skewed Brcit-Wigner distribution 



3.2.4 The 'P Distribution 

Integration of yV(cos 9, (p, <?) over cos 9 and averaging over 
(f) is equivalent to summing over the final state helicities 
(A, A'), giving 

1 



Wi<P) = — { 1 - ecos2<?>Tr(ri) + y^2e{l + e) cos<?Tr(r'^) 
27r 

(9) 



-Pb V2e(l-e) sin <P Tr (r^ ) } . 



As a sum over the p° decay states, yV(<?) measures the 
out-of-plane dependence of the eN eN reaction, and 
the matrix elements Tr(r^), Tr(r^), and Tr(r^) can be 
identified with the out-of-plane response functions famil- 
iar from medium-energy quasi-elastic electron scattering: 
Wtt , Wlt , and Wlt' , respectively (see ) . Parity con- 
servation ensures that rf^ = r"i_i, for a = 0,1,4,5,8; 
thus Tr(r°) = 2r"]^ -|- Tqq, in agreement with the notation 
of ref. If SCHC applies, then W{<P) = l/27r. 



dN 

dM^rT, 



[M, 



2 

TTTT 



(14) 



The mass-dependent width rp{MT^„) and skewing param- 
eter ns = 3.18 were taken from ref. ps[ |. 

The events were generated from an exponential dis- 
tribution in t' , and weighted to account for the coherent 
contribution at low t': 



dN _ CTA 



oc —bp^e 
(Tn 



bAt' 



(15) 



The ratio of the total coherent cross section from -^He to 
the incoherent one from nucleons was set to cta/cn = 0.48, 
and the diffractive slope for '^He (the nucleon) was taken 
as 25 (6) GeV^2 gee [3). The present data set was 
found to be consistent with (n2-n3) using the parameters 
as given above. While the value of the skewing parameter 
ns = 3.0±0.02 provides a better fit [Q, this hardly affects 
the angular distributions. 
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4.2 Procedure 

The analysis was done in four regions containing data 
of approximately equal statistical precision: 0.5 ... 0.95 
. . . 1.3 .. . 1.96 .. . 4 GeV^. For each bin, a sample 
of Monte Carlo events was generated from a distribution 
that is isotropic in all angles. These events were then sub- 
jected to the simulation of the experiment and normal 
data analysis procedures, accounting for the rather severe 
effects of the spectrometer acceptance. In order to impose 
the effects of the rj^y^, on the shapes of the simulated an- 



assuming SCHC by 

5.6,7,8 
Tl-l 

1,2,3,04 
^10 



' 1-1: 



= 0, 
= 0. 



gular distributions, the surviving Nmc 
were weighted according to 



Wi = >V(cos6'i,(/)^,^j;ei) (i = 1, ...,Nmc), 

with >V(cos6'i,0i,^i;ej) given by (|). The cos^^, 
and Ci refer to the generated kinematics of each Monte 
Carlo event. The weights Wi were evaluated as a function 
of the 23 r-matrix elements describing yV(cos 61, 0, (||) 
as follows: 

— The matrix elements that control the one-dimensional 
distributions yV(cos6'), >V(0), >V(<?), and yV(V'), were 
taken as free parameters: Tqq, r^_x^ '"i-ij ''i-ii Tr (r^), 
Tr(r5) and Tr(r8). 

— An eighth free parameter was taken to be the phase 8 
between the longitudinal and transverse Wick ampli- 
tudes Too a-iid Til — the only ones remaining under the 
assumption of SCHC and natural parity exchange |^ . 
This parameter was used to constrain the following 
matrix elements: 



The 8 free parameters plus a single arbitrary global 
normalization factor were simultaneously adjusted to min- 
imize the total of the deviations between the data and 
the Monte Carlo model in their reconstructed yield dis- 
tributions dN/dcosO, dN/d(t), dN/d^, and dN/dip. The 
1.5 X 10* events t^2 degree of freedom ranges from 29/42 to 48/42 for 
the four bins. The optimized parameters and their un- 
certainties are taken to be the measured values at each 
point. The quoted statistical uncertainties include the 
effects of the significant correlations between parameters. 



cos (5 



2 V 2 (l-f ei?) ^ 

1 pR sin 6 

{1 + eR)' 



(16) 
(17) 



4.3 Systematic Uncertainties 

Systematic uncertainties due to assumptions used in the 
analysis were estimated by performing the analysis under 
different yet tenable assumptions and comparing results. 
The following issues were investigated: 

— Geometric acceptance: The systematic uncertainty due 
to alignment errors was estimated by varying the align- 
ment in the GEANT model within known tolerances 
and performing the analysis procedure on Monte Carlo 
data. This contribution tends to dominate at small 
for all matrix elements except Tr(r^) and Tr(r^). 

— Assumptions in the VMD-based diffractive generator: 

and W dependence ((p^) and (p^)), relative con- 
tributions ((Ta/cn) from coherent and incoherent scat- 
tering and the diffractive slope parameters &n and 6a 
(^5|) were all varied within uncertainties. 

— Background correction: The systematic uncertainty at- 
tributable to the background correction was estimated 
by varying the AE cut in the analysis procedure. 



The Tin ' ' terms, due to longitudinal-transverse in- 



with R calculated from r\ 

5,6,7,8 

terference, cause a correlation between cos 9 and ^ and 
are not observable in one-dimensional angular distri- 
butions; however, the correlation between cos9 and 
in the detector acceptance function caused by the in- 
trusion of the magnet shielding plate requires that all 
matrix elements must be fitted simultaneously. 
The matrix elements describing the azimuthal produc- 
tion angle (<?) distribution were assumed to have the 
helicity structure 



and e according to (pT[). — Relative contributions of Tqq'^ and r[f''^ to the trace 



1.5,8 




and (|19|)) measured by the yV(«?) dis- 



1,5,8 
'^00 

1,5.8 



Tr(ri^^^«) rgo^ 



Trir'''^') -(1- 



'00 J J 



(18) 
(19) 



which assigns to these amplitudes the same ratio of 
longitudinal-to-transverse strength as observed in the 
distribution averaged over 

The remaining matrix elements are not observable in 
the one-dimensional distributions and were constrained 



Tr(ri^5,8) 
tribution (K 

— Possible violations of SCHC: Previous measurements 
^ covering a wide energy range have yielded non-zero 
values for tJq, which can be attributed to interference 
between a small spin-flip amplitude and a non-spin-flip 
amphtude |Tflip|/|rno-fiip| ~ 0.14 (see Appendix 1 of 
Ref. §). 

— Radiative effects: The spectrometer material in the 
path of the scattered positrons leaving the target is 
about 5% of a radiation length. This happens to be 
an amount that will produce an effect similar to that 
of internal radiation, based on the equivalent radiator 
approximation | pl[ |. Hence the systematic uncertainty 
attributable to uncorrected internal radiative effects 
was estimated by comparing fit results using Monte 
Carlo data generated with and without external ra- 
diative effects. 

The uncertainty contributions thus estimated were added 
in quadrature. 
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Fig. 7. Best-fit values for tlie pliase difference 5 between tlie 
longitudinal and transverse amplitudes (statistical uncertain- 
ties only). 

5 Results and Discussion 

The best-fit values of the phase angle S shown in Fig. ^ 
are consistent with zero, as previously observed for W > 
3 GeV [|| , and as expected when both longitudinal and 
transverse production are difFractive. 

Previous data have typically been presented as one- 
dimensional distributions with corrections for instrumen- 
tal effects applied. In keeping with this tradition, an ac- 
ceptance correction was trivially obtained from the op- 
timized matrix elements by comparing the Monte Carlo 
yields {dN/dcosO, dNjd^, dN/d^, and dN/d^p) with the 
analytic expressions (|2"|7|) for the angular distributions 
(>V(cos6'), W{(j)), Wi'P), and W(i/'), respectively). The 
corresponding measured angular distributions after ap- 
plying this correction are shown for four values of in 
Figs.|,|l|,|l| and |. These fi gures should be interpreted 
with caution, as each one alone does not adequately rep- 
resent the global fit involving strong correlations between 
the resulting matrix elements. Systematic and statistical 
uncertainties combined in the error bands apply in the 
same way to both the curves representing the results of 
the global fit and — through the acceptance correction — 
the data. 



5.1 Matrix Elements Allowed by SCHC 

As observed in previous measurements, the yV(cos6') dis- 
tribution in Fig. ^ shows a dramatic shift as increases, 
from a predominantly sin^ 6 shape to a predominantly 
cos^ 9 shape, indicating a transition from transverse to 



longitudinal production (see (|^)). The structure of yV{ip) 
visible in Fig. ^ is entirely due to transverse mesons 
{r\_i in (0)), so that as the longitudinal component in- 
creases, the 'damping' of the yV'(?/') distribution indicates 
a complementary decrease in rl_^ because of the dilution 
from longitudinal mesons. The measured matrix ele- 
ments and are plotted versus in Fig. |l^ and 
are listed in Table 0. Comparison with data of other exper- 
iments in Fig. |l^ is complicated by their different values 



of the polarization parameter e. For that reason, this issue 
will be revisited when discussing the values of R derived 
from the data in Fig. 

If SCHC is valid and the j*N p°N reaction pro- 
ceeds via the t-channel exchange of a particle with natural 
parity (J^ — 0+, 1^,...), then the tensor alignment of the 
mesons in the scattering plane is related to the transverse 
fraction of mesons by 



2rl_, = (1 



„04\ 



(20) 



Thus the decay pions from transverse p*^ mesons are pref- 
erentially emitted in the scattering plane, while those emit- 
ted from longitudinal mesons show no ip dependence. 
A comparison of the last two columns of Table ^ reveals 
that the measured values of Tqq and r\_i are in excellent 



agreement with (^0|), with a total oi 2.7 for 4 degrees 
of freedom, accounting for correlations. This demonstrates 
that the reaction is dominated by helicity-conserving am- 
plitudes with natural parity exchange. 



5.2 Matrix Elements forbidden by SCHC 

The distributions in the azimuthal decay and production 
angles and <P are shown in Figs. ^ and ^ the ma- 
trix elements describing the distributions can be found in 
Tabl es 3 and H respectively, and are plotted versus in 
Figs. O^and M. Where other data exist, they are consis- 
tent. The present data for Tr(r^) and Tr(r^) interpolate in 
between the photo- and lepto-production data from the 
collider experiments [10, but appear consistent in spite 
of the large difference in W. No previous data for Tr(r^) 



exist. The extraction of 3(rf_]^) and Tr(r°) depends on 
the use of a polarized lepton beam. Even though a single 
beam helicity was used, the effect of the term involving 
rf_i is not visible, and only a cos 2(j) variation can be seen 
in yV(0). This dependence is clearly visible in Fig. ^ at 
the smallest value of Q^, although this variation has only 
marginal significance in comparison with the total uncer- 
tainties. Within the systematic uncertainty of the mea- 
surement, the present results are in fair agreement with 
the SCHC prediction that r°t^ = 3(rf_i) = Tr(ri) = 
Tr(r^) = Tr(r^) = 0, though small negative values for r"!]^ 
and Tr(r^) (for the latter 1.8 standard deviations at the 
smallest Q^, and about 2.5 a averaged over Q^) could both 
be attributed to the interference of a small double-spin- 
flip amplitude with a non-spin-flip amplitude. Similarly, 
the non-zero value of Tr(r^) — 3.5 standard deviations 
averaged over Q^, ignoring the systematic uncertainties 
which only in this case are dominated by radiative effects 
— could be caused by interference of a spin-flip ampli- 
tude with a helicity conserviiig amplitude (see Appendix 
A of and Appendix 1 of [||). The present data are in 
agreement with the significantly positive values for Tr(r^) 
from the HERA coUiding-beam experiments at large W, 
which reflect the similar SCHC-violating behaviour re- 
ported for their values of Tqq ||i^,|11|. However, radiative 
effects are not included in their systematic uncertainties. 
This violation by Tqq was not observed in measurements 
with apparently similar precision at low 2.3 GeV) M . 
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Fig. 8. The acceptance-corrected an- 
gular distributions W(cos 6) in four re- 
gions of Q^, with statistical error bars. 
The solid curves are from evalu- 
ated with the best-fit matrix elements. 
The shaded regions indicate the uncer- 
tainty of the curves (or of the accep- 
tance correction to the points) arising 
from the total uncertainty of rfJo (see 
text). The average values of are 
shown in GeV^. 



Fig. 9. The acceptance-corrected an- 
gular distributions yV{ip) in four re- 
gions of Q^, with statistical error bars. 
The solid curves are from (^ evaluated 
with the best-fit matrix elements. The 
shaded regions indicate the uncertainty 
of the curves arising from the total un- 
certainty of rl_i. The average values of 
are shown in GeV'^. 
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Fig. 10. The spin density matrix 
elements rfJo and rl^i as a func- 
tion of (•), with the band rep- 
resenting the systematic uncertainty. 
There are shown for comparison pre- 
vious measurements with statistical 
uncertainties as (inner) error bars, 
from DESY § (v), EMC [| (A), 
NMC § W, E665 § (□), ZEUS | (o), 
and HI |hl[| (★). Where the systematic 
uncertainties were specified HP_l|, they 
are added in quadrature to produce the 
outer error bars. 







measured roo 


measured rl^i 


ri_i(SCHC) 
from r-oo 


0.80 


0.76 


0.147 ±0.030 ±0.073 


0.412 ±0.031 ±0.093 


0.43 ± 0.02 


1.13 


0.82 


0.288 ±0.040 ±0.041 


0.314 ±0.035 ±0.041 


0.36 ±0.02 


1.60 


0.81 


0.350 ±0.029 ±0.036 


0.302 ± 0.028 ± 0.031 


0.32 ±0.02 


2.75 


0.78 


0.432 ± 0.035 ± 0.034 


0.280 ± 0.034 ± 0.022 


0.28 ±0.02 



Table 2. Spin density matrix elements de- 
scribing the distributions in decay angle 
VV(cos6') and in polarization angle W{tp). 
The measured values of rl_i are compared 
with the predictions based on ( po| ) using 
the measured values of rgg. The first (sec- 
ond) uncertainty is statistical (systematic). 
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Fig. 11. The acceptance-corrected an- 
gular distributions VV(<?!)) in four re- 
gions of , with statistical error bars. 
The solid curves are from evaluated 
with the best-fit matrix elements. The 
shaded regions indicate the uncertainty 
of the curves arising from the total un- 
certainties in r^li and SJ(ri_i). The av- 
erage values of are shown in GeV^. 
The dashed horizontal line represents 
the value ^ predicted by SCHC. 
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Fig. 12. The spin density matrix ele- 
ments for rill and ^(ri.i) as a func- 
tion of (•), with the band rep- 
resenting the systematic uncertainty. 
There are shown for comparison previ- 
ous measurements with statistical un- 
certainties only from DESY [|l) (V), 
NMC § (o), E665 [| (□), and from 
HI (*) where the systematic uncer- 
tainties are specified to be negligible in 
this case. 



(GeV^) 




^1-1 


3(rLi) 


0.80 


0.76 


-0.107 ±0.036 ±0.073 


0.186 ±0.117 ±0.132 


1.13 


0.82 


-0.027 ±0.031 ±0.041 


0.134 ±0.101 ±0.086 


1.60 


0.81 


-0.035 ± 0.030 ± 0.036 


-0.119 ±0.085 ±0.049 


2.75 


0.78 


-0.009 ± 0.033 ± 0.034 


0.057 ±0.103 ±0.116 



Table 3. Spin density matrix elements describing the 
the azimuthal decay-angle distribution VV((?!)). The first 
(second) uncertainty is statistical (systematic). 



5.3 Ratio of Longitudinal to Transverse Cross Sections 

With the assumption of s-channel helicity conservation, 
knowledge of the virtual photon polarization allows ex- 
traction of the ratio of the longitudinal to transverse cross 
sections for diffractive production 



R = 



Oh 

(7t 



1 



,,04 

'oo 



i-r: 



04 



(21) 



without having to combine data at different beam ener- 
gies. The right side of (|2^) can be interpreted simply as 
the fraction rgg of longitudinal mesons normalized to 
the fraction e of longitudinal photons, divided by the frac- 
tion (1 — Tqq) of transverse p'^ mesons normalized to the 
unit flux of transverse photons. 

As in ( ^l|) depends only on the squares of am- 
plitudes, (^) is robust against SCHC violations. For in- 



stance, the SCHC-violating spin-flip amplitudes reported 
in [^||] to be as large as 14% of the SCHC-conserving 
non-spin-flip amplitude require only a 2% correction af- 
ter squaring them. Moreover, the double-spin-flip ampli- 
tudes T_ii and Ti_i couple transverse m = ±1 photons 
to m = =f1 mesons, and therefore do not affect the 
vaHdity of @). 

Equation ( pT| ) was used to determine R; the results 
are listed in Table |5[ and plotted in Fig. 15 together with 
previous data with relevant precision in the regime W > 
4 GeV. Gauge invariance restricts R to be zero at — 0, 
and dimensional arguments imply that the longitudinal 
cross section dominates at asymptotically large [ p5[ . 
All of the sets of existing data are consistent with a mono- 
tonic interpolation between these extremes, and show that 
both the transverse and longitudinal production mecha- 
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Fig. 13. The acceptance-corrected an- 
gular distributions W(#) in four re- 
gions of , with statistical error bars. 
The solid curves are from evalu- 
ated with the best-fit matrix elements. 
The shaded regions indicate the uncer- 
tainty of the curves arising from the to- 
tal uncertainties in Tr(r^), Tr(r^), and 
Tr(r*). The average values of Q are 
shown in GeV^. The dashed horizontal 
line represents the value ^ predicted 
by SCHC. 



Fig. 14. The matrix elements Tr(r^), 
Tr(r^) and Tr(r*), as a function of 
(•), with the band representing the 
systematic uncertainty. There are also 
shown previous measurements from 
HI |ll| (*), and also values from 
DESY [0] (V) (with statistical un- 
certainties only) and ZEUS @ (o), 
calculated from measurements of r\'{' 
and r^Q. Where the systematic uncer- 
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tainties were specified [^,y_l[, they are 
added in quadrature to produce the 
outer error bars. 







Tr(ri) 


Tr(r5) 


Tr(r») 


(GeV^) 










0.80 


0.76 


-0.332 ±0.085 ±0.100 


0.043 ± 0.034 ± 0.085 


-0.314 ±0.256 ±0.148 


1.13 


0.82 


-0.033 ± 0.075 ± 0.049 


0.054 ± 0.036 ± 0.048 


-0.221 ±0.215 ±0.291 


1.60 


0.81 


-0.174 ±0.083 ±0.033 


0.088 ± 0.032 ± 0.052 


0.314 ±0.186 ±0.206 


2.75 


0.78 


-0.114 ±0.092 ±0.035 


0.064 ±0.041 ±0.069 


-0.192 ±0.243 ±0.150 



Table 4. Spin density ma- 
trix elements describing the 
production angle distributions 
The first (second) un- 
certainty is statistical (sys- 
tematic). 



nisms are important at the energy and accessible at 
Hermes. 

The constraint R{Q^ = 0) = and the prediction by 
several models that the longitudinal and transverse cross 
sections should differ by some power of [|^,|l^,^ sug- 
gest a fit of the form 



i?(Q2) = co(W) 



(22) 



The available data now span a wide range of centre-of- 
mass energies both above and below W = 7GeV. This 
is close to the energy at which the VF-dependence of the 
production cross section makes a transition from a de- 
creasing behavior attributed to Reggeon exchange to a 
1^0.22 dependence attributed to Pomeron exchange. The 
improved precision that the present data add to the lower 
energy region allows the fit to be done with two values of 



Co |12[ : one for the data from high energy muon beams 
and collider experiments spanning 7 < < 140 GeV 
(EMC §, NMC §, E665 §, ZEUS @, and HI 
and another for the present data in the lower beam en- 
ergy range 3.8 < < 6.5 GeV, yielding 



Co 



0.32 ±0.04, 
0.48 ±0.03, 

ci 0.66 ±0.03. 



4 < (VF) < 7 GeV 
{W) > 7 GeV 



The P^r degree of freedom for the combined fit is 
20.7/20, calculated from the statistical uncertainties, plus 
the systematic uncertainties in quadrature where the lat- 
ter are available. Difficulties in assessing the systematic 
uncertainties in some of the data sets preclude a compre- 
hensive treatment of all the systematic uncertainties. 

The only previous R data |l^ in the Hermes energy 
region is consistent but insufficiently precise to influence 
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Fig. 15. The ratio R = ctl/ctt determined from tiie measured 
values of t-qo using (^l|) (•) . The errors bars are statistical only, 
while the shaded region indicates the systematic uncertainty of 
the present data. Previous measurements are shown for com- 
parison, with experiment names as in Fig. |l^ with the addition 
of Cornell for ref. . Where the systematic uncertainties were 
specified |^ 1^ , they are added in quadrature to produce the 
outer error bars. The dashed (dotted) line represents a simul- 
taneous fit according to ( |2^ ) to data with average W > 4 and 
above (below) 7 GeV, which are shown as open (filled) data 
points. The DESY and CorneU data for 2 < < 4 GeV are 
also plotted but were not included in the fit (see text). 

Table 5. i? as determined from Tqq. 



(GeV^) 




R 


0.80 
1.13 
1.60 
2.75 


0.76 
0.82 
0.81 
0.78 


n. 90+0.06+0.15 
U.^O_Q.05-0.12 
r, ^Q+0.10+0.10 
U.y:»_o.09-0.09 

n flft+o.og+o.ii 

U.UU_Q Qg_Q j^g 

n QQ+O.15+0.14 
'-'■'^"-0.13-0.13 



the fit. Inclusion of either or both of two existing data 
sets [|l|,H in the range 2 < W < 4: GeV changes the results 
of the fit by less than one standard deviation, and does 
not significantly improve the precision. (It may be noted 
that the data oif reference Q are not well-described by the 
fit function — their contribution would be similar to 
that of all of the rest of the fitted data in both energy 
regimes.) 

As observed in previous experiments, the depen- 
dence of R is in marked disagreement with the predic- 
tion of either VMD || or pQCD H that a = 1. The 
value of cq is found to increase with centre-of-mass energy 
in the range above W = A GeV where both longitudinal 
and transver se p roduction are clearly diffractivcH It has 
been shown l25| that the energy dependences of the lon- 



^ Experiments at lower energy |l|,|^ provide evidence that 
for fixed Q^, R increases also as W falls below 2 GeV. 



gitudinal and transverse couplings for gluon exchange are 
identical, so that R is independent of energy at high en- 
ergy; however, recent theoretical work ||l^ has suggested 
that at lower energy where quark exchange is dominant, 
there may be an additional class of diagrams that enhance 
the transverse couplings. The observed energy dependence 
of Co is consistent with this trend. However, it should be 
kept in mind that the lower energy fit to the present data 
set on the ■^He target could be influenced by nuclear ef- 
fects, although it has been suggested that this would in- 
crease R Q . Other matrix elements accessible with a real 
photon beam were found to be the same for coherent p° 
production on deuterium and incoherent production on 
either nucleon in deuterium . Furthermore, those data 
are consistent with other similar measurements on hydro- 
gen iQ. Finally, coherent production on such a light nu- 
cleus as ^He accounts for only about 25% of the events (see 
Fig. Hence a nuclear effect is unlikely to account for 
the observed difference from the results at higher energy.^ 



6 Summary 

The production and decay-angle distributions for the com- 
bined coherent and incoherent polarized electroproduction 
of mesons from '^He have been measured at photon vir- 
tuality 0.5 < < 4 GeV^ and photon-nucleon centre-of- 
mass energy 3.8 < T4^ < 6.5 GeV. The spin-density matrix 



elements rgg. 



„04 

"i-ij 



_i, Tr(ri), Tr(r5) and Ti (r^) 
are found to be mostly consistent with s-channel helicity 
conservation and natural parity exchange in the i-channel. 
Possible deviations from SCHC at the few sigma level were 
found for Tr{r^). The ratio R of the cross sections for lon- 
gitudinal and transverse photons was determined from the 
decay-angle distribution over a range 0.5 < < 4 GeV^ 
and at (W) = 5 GeV. The dependence of R is consis- 
tent with previous experiments, exhibiting a monotonic 
increase with Q^; however, the present data now sug- 
gest that at a given value, R increases with average 
photon-nucleon centre-of-mass energy for such energies 
above 4 GeV, assuming that nuclear effects in coherent 
production on ^He can be neglected. 
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